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Abstract 

The dynamics of one dimensional iterative maps in the regime of 
fully developed chaos is studied in detail. Motivated by the obser- 
vation of dynamical structures around the unstable fixed point we 
introduce the geometrical concept of a turning point which represents 
a local minimum or maximum of the trajectory. Following we in- 
vestigate the highly organized and structured distribution of turning 
points. The turning point dynamics is discussed and the correspond- 
ing turning point map which possesses an appealing asymptotic scal- 
ing property is investigated. Strong correlations are shown to exist 
for the turning point trajectories which contain the information of the 
fixed points as well as the stability coefficients of the dynamical sys- 
tem. For the more specialized case of symmetric maps which possess 
a symmetric density we derive universal statistical properties of the 
corresponding turning point dynamics. Using the turning point con- 
cept we finally develop a method for the analysis of (one dimensional) 
time series. 
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1 Introduction 



In the past twenty years one dimensional iterative maps have attracted much 
attention and became a very active field of research. An important reason 
for the strong interest in the dynamics generated by iterative maps is the 
fact that they exhibit a variety of typical nonlinear phenomena such as pe- 
riod multiplication, intermittency and chaotic behaviour [Grossmann and 
Thomae, 1977, Feigenbaum, 1978, Feigenbaum, 1979, Collet and Eckmann, 
1980, Crutchficid ct aL 1982, Ott, 1993, Schuster, 1994]. They provide there- 
fore a tool for the modeling and simulation of the dynamics of more compli- 
cated systems as they typically occur in physics, chemistry and biology. In 
particular the chaotic dynamics has by now been shown to be of potential 
importance in many different fields like for example fiuids [GoUub and Ben- 
son, 1980], plasmas [Sagdccv et al, 1990], circuits [Linsay, 1981] and lasers 
[Arecchi et al, 1982]. Many of the characterstic features observed in dissi- 
pative systems of higher dimensional phase space have their origin in the 
universal behaviour of one dimensional processes. 

Of particular relevance are the maps of the unit interval which possess a 
single smooth maximum. Due to their noninvertibility they show a very rich 
behaviour when a control parameter is changed continuously. A typical phe- 
nomenon is the bifurcation route to chaos and its universal scaling laws which 
have been studied in detail in the literature [Grossmann and Thomae, 1977, 
Feigenbaum, 1978, Collet and Eckmann, 1980] and are now well understood. 
This bifurcation route leads to an attractor which is, with further changing 
control parameter, followed by an inverse band bifurcation and ends up in 
the final state of fully developed chaos. It is the purpose of the present paper 
to investigate the dynamics of these maps for fully developed chaos, i.e. in 
their chaotic and ergodic state. Common quantities describing this state are 
the Ljapunov exponent A which yields the rate of exponential divergence for 
neighbouring trajectories and the invariant density p which is a measure for 
the distribution of orbits in the ergodic limit. For further characterization 
one can use the non-uniformity factor which measures the deviation of the 
rate of divergence from its average or the correlation function [Gyorgyi and 
Szepfalusy, 1984]. Our investigation is aimed at a further understanding of 
the state of fully developed chaos and goes beyond the above quantities [Di- 
akonos and Schmelcher, 1997]. This will provide us valuable insights into the 
underlying structure of the ergodic dynamics of maps. We emphasize that 
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the turning point properties and concepts developed in the present paper 
are vahd for any one dimensional smooth and single humped map for which 
the rank two image of the extremum merges into an unstable fixed point 
with a positive multiplier. We shall call this rather general class of maps 
in the following SSH maps. The only exception is the universal statistical 
behaviour derived in Sec. 5 which holds only for doubly symmetric maps, i.e. 
for symmetric maps which possess a symmetric density. 

In detail we proceed as follows. In chapter 2 we begin by illustrating 
some eyecatching features of a typical one dimensional ergodic trajectory. 
These features do not show up in the invariant density but are related to 
the distribution of the oscillations of the trajectories. We also introduce in 
chapter 2 the concept of turning points and discuss the role of turning points 
in the ergodic dynamics. A comparison of the chaotic map with a random 
map with the same invariant density clearly reveals the structures embedded 
in the fully chaotic and ergodic dynamics. In chapter 3 the turning point 
map (TPM) is introduced and analyzed. It possesses an important scaling 
behaviour which develops by joining together humps of successive iterates of 
the original map. Chapter 4 is devoted to a turning point analysis of higher 
iterates of the map. Important dynamical information is shown to be present 
in these higher TPMs. In particular higher periodic orbits, their preimages 
as well as local stability coefficients are exhibited in the TPMs. Chapter 5 
contains a derivation of the universal statistical turning point properties for 
doubly symmetric maps. In chapter 6 we use our turning point concept in 
order to develop a method for the analysis of one dimensional time series. 
Chapter 7 concludes with a summary. 

2 The dynamics and the distribution of turn- 
ing points 

In the ergodic case of fully developed chaotic dynamics the measure of the 
chaotic trajectories is equal to one. Unstable periodic orbits are dense in this 
chaotic 'sea' but they are of measure zero. Let us begin our investigation by 
considering a typical chaotic trajectory. For reasons of illustration we hereby 
choose the logistic map in the ergodic limit x{n + 1) = 4x(n)(l — x{n)). We 
emphasize however that our main observations and results of investigations 
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are valid for any SSH map. Figure 1 shows such an ergodic trajectory for 
1.5 X 10^ iterations together with its invariant density for 10^ iterations. The 
invariant density of the logistic map is given by pi, = — , Both its 

■K^ x{l — x) 

smoothness in the open interval ]0, 1[ as well as its singular behaviour at 
X = 0, 1 can clearly be seen in the density of the trajectory given in Fig.l. 

A closer look at the trajectory itself in Fig.l gives us the impression that 
there exists a dynamical structure in the neighbourhood of = | which 
is the position of the unstable fixed point {xp refers in the following always 
to the fixed point different from zero of the map). This structure possesses 
no counterpart in the corresponding smooth invariant density Pl{x). Figure 
2 shows the trajectory for fewer iterations which reveals that the above- 
mentioned structure in the neighbourhood of the fixed point might have 
its origin in the nonuniform distribution of oscillations in the interval [0, 1]: 
we observe an enhanced probability for finding small amplitude oscillations 
in the vicinity of the fixed points. Before investigating these properties in 
more detail let us introduce the following notions: a turning point x(n) of a 
trajectory (abbreviated by TP) is a point which represents a local maximum 
or minimum of the trajectory (see Eq.(l) below for a quantitative criterion 
for the turning point). We remark that the name 'turning point' should 
not be confused with the same name that often refers to the local extrema 
of a map /. The center of oscillations (abbreviated by CO) are defined to 
be the midpoints between subsequent local maxima and minima (and vice 
versa) and the abbreviation AO means the amplitude of the corresponding 
oscillations. 

In the following we analyze the dynamical properties of the turning points. 
This will lead us to a qualitative as well as quantitative understanding of 
the mechanisms which create, among others, the corresponding structures 
(see below) in the relevant densities. We begin by deriving a necessary and 
sufficient condition so that a point x{n) of an ergodic trajectory is a turning 
point. Our starting point is the geometric condition 

{x{n + 1) — x{n)) X {x{n) — x{n — 1)) < (1) 

for a turning point x{n) where x{n + 1) = f{x{n)) and / is a single humped 
and smooth map. Xi is the inverse image of the fixed point xp on the left 
branch [0, Xm] of the map / where x^ is the maximum of the map /. We 
will show in the following that given a map / and its fixed point xp the 
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condition for x{n) to become a turning point is x{n — 1) > Xi. Since the 
map / is monotonically increasing on [0,Xr„] and decreasing on [x^, 1] we 
have to distinguish three cases. For x{n — 1) < Xi we obtain the inequahty 
x(n — 1) < x{n) < x{n + 1) and therefore x{n) is not a turning point. In 
contrast to this we obtain for x(n — 1) G [ ] the inequahties x{n — 1) < 
x{n) and x{n) > x{n + 1) and x{n) is a turning point. The third case is 
x{n — l) e [xi?, 1] which yields a; (n — 1) > x{n) andx(n + l) > x(n), i.e. x{n) 
is a turning point. In total we therefore arrive at the condition x{n — l)>xi 
for x{n) to become a turning point. Otherwise, i.e. for x{n — 1) < xi, x{n) 
is a non turning point. For the special case of the logistic map this means 
x{n — 1) has to be larger than j in order for x{n) to become a turning point. 

With this result we are now able to understand and discuss the occurence 
of turning points, i.e. we gain insights into the dynamics of turning points, 
(see in particular also Fig. 2 for an illustration of the turning point dynamics). 
In the following we distinguish between different phases of the motion created 
by the iterations of the map. Successive iterations of the map which are 
located in the interval [0,xi] belong to a stretching phase of the iterative 
motion. Apart from the very first point (see below) in the interval [0,xi] 
the points x{n) of this phase obey x{n — 1) < xi and are therefore non 
turning points. Once we leave the interval [0, Xi] the first point x > Xi obeys 
X G and is also a non turning point. The next series of iterated 

points obey alternately x > xp > xi and xi < x < xp and we therefore 
obtain an oscillating stretching phase of motion around the fixed point xp 
which consists exclusively of turning points. These stretching phases end in 
case we obtain x > x* where x* is the inverse image of Xi on the right branch 
[xm, 1] of the map / (for the logistic map we have x* = | + |(v^ — 1)). The 
next iterated point x is then folded back into [0, Xi] and represents a turning 
point. Subsequently the next stretching phase of non turning points in the 
interval [0, Xi] starts etc. Both the stretching phase as well as the oscillatory 
phase of motion can nicely be seen in Figure 2. 

In order to gain further insights into relevant turning point properties we 
have illustrated in Fig. 3 the density prp of the turning points for a typical 
chaotic trajectory of the logistic map. First of all we realize that prp exhibits 
a step like structure at the position x^ of the unstable fixed point. We observe 
an enhanced probability of finding turning points in the interval [xp, 1]. The 
singular behaviour of the turning point density at x = and x = 1 is, apart 
from a constant factor, the same as that for the density and characteristic 
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for the underlying dynamical law. 

It is possible to understand the origin of the step-like structure of the 
turning point density ptp in Fig. 3 analytically. Using the above condition 
x{n — 1) > Xi for x{n) to become a turning point we obtain the density of 
the preimages of the turning points: M Pl{x)Q{x — Xi). Here is the step 
function G(x) = {1 for a; > 0; for a; < 0} and M the corresponding 
normalization constant. Mapping this density forward with / we get the 
expression ^pl{x){\-\-Q{x — xp)) for the density of the turning points. This 
short derivation clearly shows that the step like structure of the turning point 
density at the position of the fixed point is a common property for all SSH 
maps. Looking at the turning point density therefore immediately reveals 
the positions of the fixed points different from zero for a given SSH map. 

In order to get an impression of the dynamical information contained in 
Fig. 3 we have to compare the turning point density pxp of the logistic map 
with the density of turning points Prtp which is generated by the weighted 
random map. (The weighted random map is chosen to possess the same 
density as the chaotic map, i.e. pr = pi)- 

Indeed the density pbtp can be derived in closed form via the following 
considerations. Given the tent map g{y) — {2y for ?/ < |; (2 — 2?/) ify>^} 
it is well-known that the conjugacy between g and the logistic map / is given 
by the homeomorphism h(y) = kz^Slilll^ [ q ^ g — ^ j^^ ^^le invariant 
density of g is equal to one on the whole unit interval. If y{n) = y is a. point 
of the interval, the probability that it is a turning point is y'^ + {1 — y^. 
The normalized random turning point distribution ajixp belonging to the 
random map with unit density, i.e. 0"^ = 1, is therefore given by cfrtp = 
^[y^+{l—y)^]. The enhancement of the turning point density apTP compared 
to (jR at and 1 are geometric boundary effects. Since the conjugacy h 
preserves the turning points and transports the invariant density for g to 
the invariant density of / we can treat the application of /i as a coordinate 
change of the unit interval. As a result the weighted random turning point 
density Prtp takes on the following appearance 




CTRTpiv) 



6[{arcsin^/xy + {arccos^/xy] 






Prtp is symmetric with respect to | and smoothly monotonically decreasing 
and increasing in ]0, |] and [|, 1[, respectively. It possesses power singularities 
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at X = and x = 1 which are apart from a constant factor the same as for 

PR- 

To investigate this further we have illustrated in Fig.4a the density for 
the centers of the oscillations Plco and in Fig. 4b the density for the ampli- 
tudes of the oscillations Plao for the trajectory of the logistic map shown 
in Figs. 1/2. First of all we observe that the density pico is nonzero only 
in a subset of the unit interval, namely for approximately x G [0.38,0.78]. 
The density for the centers of the oscillations shows three main peaks. The 
largest peak is located at |. Two other major peaks are located at ap- 
proximately 0.58 and 0.78, respectively. The centers of the oscillations are 
therefore not smoothly distributed but exhibit certain dynamically prefered 
values. In addition the density Plco shows a number of small dips and step 
like structures. The density for the amplitudes of the oscillations shown in 
Fig. 4b exhibits a progression of four sharp peaks which arc preceded by step 
like and strongly ascending broad structures. The positions of the peaks 
are approximately 0.56,0.86,0.96 and 1.0, respectively, which represent the 
dynamically prefered amplitudes of the oscillations. Both the distribution 
Plco of the centers as well as the density Plao for the amplitudes of the os- 
cillations of the trajectories contain therefore dynamical information which is 
in particular responsible for the above-observed structures in the trajectory. 
We emphasize that looking at both Plco,Plao with increasing resolution 
reveals an increasing number of structures and detailed properties. Part of 
these structures can be derived analytically by the following considerations. 

Consider the points xo — 1 and Xk — (/|[o,1])~'^(^f) for A; > which are 
the preimages of the fixed point on the left branch of the map /. If the point 
x(n) = X is a turning point and x G {xk, x^-i] then the next turning point is 
x{n + k) = f^{x). Set 

^.(^) = ^^^4^ 
which defines the center of an oscillation. Then 

oo 

Plco = Y.{fk)*{pTPX{xk,Xk-A) (4) 
k=i 

where xa denotes the characteristic function of A (1 on A and outside A). 
The Perron-Frobenius operator {ipk}* transforms the densities in the usual 



7 



way 



E 



pjy) 
\v'k{y)\ 



(5) 



In particular, the support of Plco is the union of the supports of the densities 
{fk)*{PTPXixk,xk_i]), that is the union of the sets {'^k){{xk, Xk-i])- We have 
therefore for the logistic map 



ipi{x) 



5x — ix^ 



(6) 



and {xq, Xi] = (^,1]. The critical point, i.e. maximum, of (pi is |. In addition 
V^i(i) = i "PiiD = 1' and <^i(l) = i. Therefore v?i((xi, Xq]) = [|,§]. 



(7) 



If A; > 1 and x e (xk, Xk-i] then f'^ ^{x) e {xi, xp] and ■'^ > x so 

<^fcW = 2 2 " '^^^■^ ^"^^^ - 32 

On the other hand, /'^(x) e [a:;^, 1], so 

/*^(a;) +x fix) xf 3 

Therefore the support of pLco is contained in [|, ||]. Moreover, for /c > 1 
there are points ak,bk G {xk,Xk-i] with f^{ak) — xp and f^ipk) = 1. We 
have 



(8) 



C ipk{(xk,Xk-i]) 



(9) 



and 



, (a;fc_i+a:F) b^+l 1 
— 9. ' 9. ^ 9.1 



SO 



2^2 
Xk-1 +Xf I 



C (^fe((a;ik,a:fe-i]) 



(10) 



Therefore the support of pico is equal to [|, gf]- 

Apart from the support of the quantity Plco also its singularity structure 
can be derived. Plco possesses singularities whenever pxp has a singularity or 
ip'^ is zero. There are two singularities of Ptp'- at and 1. The singularity at 1 
gets carried by to i, so Plco has a singularity at |. The singularity at 
gets "washed out" , since as A; — > oo, the derivative of pk is of order 2^, whereas 
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Ptp on (xfc,Xfc_i] is only of order 2~ . However, the zeros of Lp'j^ produce 
additional singularities of pLco (one singularity for every k). In particular, 
'^i(l) — 0' get a singularity of Plco at ||. For k — 2we have <^2(t) = 

for t ^ 0.163 (t is a solution to 256t^ - SSU"^ + 160t - 17 = 0) and wc get a 
singularity of pLco ^2{t) ~ 0.577. To summarize, the distribution for the 
center of oscillations pico possesses a highly nontrivial structure on a finite 
subset of the unit interval and shows in particular infinitely many singularities 
which reflect the fact that the inflnitesimal turning point neighbourhood of 
the origin is mapped to arbitrary detail to the turning points in the finite 
neighbourhood of the fixed point xp (see below in the context of the turning 
point map). 

To make the structures contained in the ergodic chaotic dynamics even 
clearer we have illustrated in Fig. 5a and 5b the corresponding densities Prco 
and Prao for a weighted random map which possesses the same density pu — 
Pl as the logistic map. The density for the centers of the oscillations pRco 
of the random map exhibits a single broad peak which is centered around 
I and decays monotonically towards the boundaries of the unit interval. 
The most probable value for the position of the center of the oscillations 
is therefore | and the probability for values different from | decreases with 
increasing distance of the corresponding value from |. In comparison with 
the density pico (sec Fig. 4a) of the logistic map whose origin is a dynamical 
law the density pRco belonging to the random map shows a very simple and 
smooth structure. Figure 5b shows the density Prao for the amplitudes of 
the oscillations for the same weighted random map. Prao is monotonically 
increasing on the unit interval with a strongly peaked maximum at x = 
1. Within the statistical accuracy of Fig. 5b Prao is a smooth distribution 
without any structures. This clearly reveals that the fully chaotic system 
contains dynamical information according to the distributions Plco and Plao 
which is not present in the corresponding distributions Prco and Prao the 
random map. 

To complete our picture of the oscillatory motion of the random map on 
the one hand and of a ID dynamical law on the other hand we have illustrated 
in Fig. 6 the occupied coordinate space for the amplitudes of the oscillations 
as a function of the center of the oscillations for both the weighted random 
map xrao{xrco) in Fig. 6a as well as the logistic map xlao{^lco) in Fig.6b. 
For the case of the random map, xrao{xrco) AUs the two-dimensional area 
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of an isosceles triangle (see Fig.6a). In contrast to this the corresponding 
quantity xlao{xlco) for the logistic map is a one dimensional curve (see 
Fig. 6b) which consists, apart from a cross diagonal hne, of a series of skew 
and increasingly peaked humps. With increasing absolute values of these 
humps their widths tend to zero, i.e. the series of humps approaches an 
accumulation line which is the limiting line of the upper left part of the 
curve shown in Fig. 6b. Again it is obvious that there exists a major difference 
between the 'motion' of the random map which is distributed according to 
its weight and the highly organized motion of the logistic map which shows 
according to Fig. 6b a number of interesting structures. 

We emphasize that the above-observed properties and behaviour of the 
oscillatory motion of maps are by no means restricted to our specific exam- 
ple, namely the logistic map: all SSH maps show a very similar behaviour. 
This includes the special family of symmetric beta maps which has been 
introduced recently [Diakonos and Schmclchcr, 1996] and possesses mem- 
bers with a wide range for the order of the maximum and for the singular 
behaviour of the densities at the positions x = 0, 1. 

Let us now investigate how subsequent turning points are mapped unto 
each other. 

3 The turning point map and its scaling prop- 
erties 

Having discussed the occurence of turning points in the dynamics of the map 
/ we introduce now the turning point map (TPM), i.e. the map of a turning 
point onto its subsequent turning point. The TPM will be a key quanitity for 
our understanding of the structures and organization in the ergodic dynamics 
of iterative maps. Figure 7 shows the TPM belonging to the logistic map on a 
logarithmic scale. For reasons of illustration we have chosen the logistic map 
as a specific example. However, as already emphasized our arguments are 
valid for any SSH map. The monotonous part of the TPM for xp < x{n) < 1 
has its origin in the fact that the turning points in the interval [x^, 1] are 
mapped via the TPM onto the interval [0,^^;']. The first hump in Fig. 7 
represents the image of the interval [xijX^] with respect to the TPM. Since 
turning points in the interval [xi, 1] are mapped directly onto turning points 



10 



by the map / itself the TPM in [xi, 1] simply reflects the shape of / in this 
interval. 

The eyecatching feature of the TPM for a; < Xi is the extreme similarity 
of its subsequent humps with decreasing values of x{n) on a logarithmic 
scale. Approaching smaller and smaller values for x{n) Fig. 7 suggests an 
asymptotically exact scaling law for subsequent humps. The TPM possesses 
an accumulation point of its humps at the origin and has dimension 1. For 
values X < xi the TPM shows us how turning points in the neighbourhood 
of the origin C/e(0) (which have their preimagc in the neighbourhood Us{l) 
of 1) are mapped onto the turning points in the finite interval [xi?, 1]. The 
TPM docs this according to Fig. 7 to arbitrary detail and in a self repeating 
form if we approach the origin. 

As indicated in Fig.7 the vertical positions on the axis x{n + 1) of the 
turning points {xj} = (xi, X2, X3, X4, ...) of the TPM are all the same and 
given by the fixed point xp. The horizontal positions of these turning points 
on the axis labeled x{n) are the positions of successive preimages of the fixed 
point Xp on the left branch of the map /. The abscissae of xi is therefore the 
preimage f'^^xp) of the fixed point xp in the interval [0,0;^], the abscissae 
of X2 is the preimage f~^{xi) of xi in the interval [0,Xm], etc... This series 
of preimages obviously converges towards the unstable fixed point and the 
ratio of subsequent preimages converges towards the derivative of the map / 
at the origin, i.e. we have 



For the special case of the logistic map we havexi = |(l-|y2 + y2 + ....V2 + x/3) 

with i — 1 nested roots. 

But not 'only' the turning points of the TPM {xi} scale asymptotically 
with a factor n but in particular also successive humps defined in the intervals 
Aj = [xj+i , Xi\ converge asymptotically with the same ratio against a final 
shape. A numerical as well as analytical closer investigation of the i-th hump 
of the TPM in the interval Aj reveals that it is identical to the top of the first 
hump of the i-th iterate of the map /. More precisely: the i-th hump of 
the TPM is identical to the intersection x{n > xp of /' which is closest 
to the unstable fixed point 0. Putting everything together the TPM is given 
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by 

ftp = r^\x) VxgA,, Vz>0 (12) 

with Ao = [xi,l]. Approaching the origin the TPM can therefore be con- 
sidered as a sequence of the top of the humps of higher and higher iterated 
maps. We remark that even the position of the humps of the TPM is the 
same as for the corresponding iterates /*. 

In addition the series of humps of the TPM converges, apart from the 
scahng, against an asymptotic final shape. This asymptotic form can be 
obtained by using the scahng property of the humps f^~^^{x) = f"'{Snx) 
which are now defined on the interval < x < Sn+i with being the first 
zero, i.e. the zero closest to the origin, of the iterate 5„ = is the 
ratio of the first zeros of successive iterates of the map /. Introducing a new 
rescaled variable z — which is defined on the whole unit interval the 
renormalized asymptotic functional form of the complete first hump is given 

by 

M^)=]imrM (13) 

This equation provides a method for the systematic approximation of the 
asymptotic form h{^z) by calculating successive iterates of the map as well as 
their first zeros. The convergence properties of this procedure turn out to be 
excellent. 

Having understood much of the dynamics of the turning points and in 
particular the structure of the turning point map we turn now to the ques- 
tion of the correlation of the dynamics of the turning points. It is well- 
known [Gyorgyi and Szepfalusy, 1984] that the strongest decay of the corre- 
lation function occurs for the ergodic limit of the doubly symmetric maps. 
They possess 5-correlated processes, i.e. their subsequent steps are com- 
pletely uncorrelated. Using the common definition of the correlation func- 
tion C/(n) = ((a; - (a;)/) (/"(a;) - (x)/))/ this means Cf{n) = C/(0)(5„o- Fig.8 
shows now the absolute value of the correlation function of the dynamics of 
the turning points in comparison with the corresponding correlation function 
of a typical ergodic trajectory for the case of the logistic map. As expected 
the correlation function of the chaotic trajectory decays completely with the 
first time step. (For a non doubly symmetric map the generic decay behaviour 
can be much 'slower'). In contrast to this, however, the correlation function 
of the turning points shows only within the first two steps a minor decay and 
stays subsequently constant ! This reveals that the turning points obey a 
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highly correlated dynamics with infinite correlation length thereby suggest- 
ing the interpretation of the turning points as a constituent (phase) of the 
dynamics which is at a critical point of a phase transition [Beck and Schlogl, 
1993]. 

Let us briefiy summarize and comment on our previous results. Our in- 
vestigation and discussion of the turning point dynamics, the turning point 
map as well as the correlation function clearly show that a turning point 
analysis extracts valuable dynamical information from a given ergodic and 
fully chaotic system. The TPM contains details of the iterates /* for all i 
and maps in particular arbitrarily small intervals in the neighbourhood of 
the origin onto the finite interval [zi?, 1]. It obeys an asymptotic scaling law 
and converges to a repeating functional form which can be derived accord- 
ing to the above-given construction scheme. The fixed points as well as its 
preimages ,r^' (see in particular also chapter 4 which provides an investiga- 
tion of the turning point maps for higher iterates of the map /) play hereby 
a particularly important role. The turning point dynamics shows an infinite 
correlation length of the corresponding correlation function, i.e. the turning 
point analysis extracts the correlated part from the original trajectory. The 
origin of this strong correlation is the oscillating dynamics around the fixed 
point xp. 

In extension of the construction scheme for the TPM it is possible to 
define so-called renormalized limit functions which are the sequence of suit- 
ably scaled humps of the infinitely high iterates of the map /. These hmit 
functions include not only the first hump (see Eq.(13)) but all oscillations of 
the infinitely high iterates lim /" of the original map /. To achieve this we 

n— >oo 

define 

fki^) = n{rt^ - rtl,)z + rt\) for . G [0, 1] (14) 

where r^'^'' is the k-th zero of the corresponding iterate Equation (14) 
scales neighbouring zeros which converge to the same value with increasing 
iterations to the unit interval. The desired renormalized family of limit func- 
tions g^'^\z) which describe the humps between the zeros r^}^ and r^."^ for 
n ^ oo, i.e. of the infinitely high iterate of /, is then given by the limit 
process 

^«(^) = Jim r((ri") - ra)z + r^) for z e [0, 1] (15) 
Let us illustrate the above procedure for the example of the logistic map in 
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the ergodic limit. The n — th iteration x{n) of a point a;(0) can for this case 
be given analytically [Pincherle, 1920, Kadanoff, 1983] 



x{n) = :j [1 - cos {2''arccos (1 - 2x(0)))] 



(16) 



The zeros r["^^ take on the following appearance 



1 — cos 



kn 



(17) 



Applying the renormalization and limiting procedure according to Eq.(15) 
the following result for the limit function g^^^ (z) of the logistic map can be 
derived 



9l \z) = - 



1 - cos 27rV(A; - 1) + {2k -l)z 



for k= 1,2,... and z e [0,1] 
(18) 

The next chapter is devoted to a turning point analysis of doubly iterated 
maps p and will give us a general idea of the information contained in 
arbitrarily high iterates 



4 Dynamical information of the turning point 
maps for higher iterates of the map 

Let us begin our investigation of the TPMs belonging to higher iterates 
by deriving the general turning point criterion. According to Eq.(l) turning 
points occur if 

A{x{n - 1)) = {F{F{x{n - 1))) - F{x{n - 1))) x {F{x{n -1)) - x{n-l)) < 

(19) 

where F = can now be any iterate of the original map /. In order to 
determine the transition from turning points to non-turning points or vice 
versa in the criterion Eq.(19) we have to look for those zeros x{n—l) of A for 
which A changes from positive to negative values or vice versa, i.e. the zeros 
with ^(x(n— 1)) 7^0. If a:(n—l) is a fixed point we have both ^(a;(n—l)) = 
and ^(x(n — 1)) = and x{n — 1) is not a point which separates turning 
points from non turning points with respect to x{n) = F{x{n—1)). However, 
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if x{n — 1) is the prcimagc of a fixed point and consequently x{n) a fixed 
point we obtain A{x{^a — 1)) = and in particular 

^{x(n - 1)) = (xin - 1) - xin)) x F\x{n - 1)) x (1 - F'ixin)) ^ (20) 
ox 

since \F'[x[n))\ > 1 (we have only unstable fixed points) and we consider only 
those preimages of the fixed points which possess a derivative F\x{n — 1)) ^ 
0. This means that the preimages of the fixed points determine the transition 
points from turning points to non turning points or vice versa. We therefore 
obtain the following turning point criterion. The preimages {p-i} of the fixed 
points define a unique partitioning of the unit interval [0, 1] into subintervals 
h — bi)Pi+i]- -l-^ order for x{n) to become a turning point x{n — 1) has to 
obey the necessary and sufficient condition x{n — 1) € U^^Q/2i+i. 

Let us illustrate this for the example of the twice iterated logistic map 
F = P- P possesses four fixed points, namely {xpi} = {0, ^^^^ |} 
with alternating positive and negative derivatives P'{xpi), respectively. The 
preimages of the fixed points which obey the above condition F'{pi) ^ are 




1 ll5 + V5 
= (2-2V^^ 



P3 = 
P5 = 




2' 2 



P9 = (o + oV^^) (21) 



'2 2 



and in addition we have for the first and last point po = and piQ = 1, 
respectively. x{n — 1) has therefore to be located in one of the intervals 
/i, /a, /g, /y or /g if we want x{n) to become a turning point. 

Before beginning with the discussion of TPM's of higher iterates /" some 
remarks are appropriate. There exist two different classes of fixed points 
for our map The first class contains the fixed points with a positive 
derivative, i.e. P' {xpi) > 1. The local dynamics around these fixed points 
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can be characterized as an unstable stretching behaviour. The dynamics in 
the neighbourhood of these points contains no turning points. The second 
class consists of the fixed points with a negative derivative, i.e. f'^'{xFi) < 
— 1. The local dynamics around these fixed points can be characterized as 
an alternating stretching behaviour which means that successive points are 
on opposite sides of the corresponding fixed point. Each point of this local 
dynamics is a turning point. Having these facts in mind we turn now to our 
discussion of the TPM of doubly iterated maps. 

Figure 9 shows the TPM of the doubly iterated logistic map on a loga- 
rithmic scale. First of all we observe that there are a number of disconnected 
objects each of which exhibits a detailed and relatively complicated structure. 
For x{n) > we observe the objects labeled by (6), (c) and (d) which ap- 
pear only once. In contrast to this we observe for x{n) < yz with decreasing 
values of x{n) a repeating sequence of structures. The absolute unit of these 
repeating structures consists of two non-overlaping objects which are labeled 
in Fig. 9 by (a) and (e), respectively. Since the repetition of these structures 
on the logarithmic scale becomes asymptotically exactly periodic we obtain 
an asymptotic scaling behaviour with respect to the above-mentioned abso- 
lute unit consisting of the asymptotic form of the objects (a) and (e). This 
scaling behaviour is similar to the one observed in the previous chapter for 
the TPM of the map / itself. The scaling factor is the derivative of the map 
P at the origin namely p (0), i.e. the stability of the fixed point 0. Be- 
fore entering the discussion of each of the objects (a), (6), (c), [d] separately, 
we want to extract some important properties from the TPM of the doubly 
iterated map p which can easily be extended to arbitrary iterates /*. As 
indicated in Fig.9 the vertical positions of the lower and upper horizontal 
edges of the objects (a), (6), (c), ((i), (e) etc. are given by the positions of 
the fixed points XFiiXFi^xpz^xpi and the point 1. The valuable information 
of the positions of the fixed points can therefore easily be extracted from 
the corresponding TPM. Higher periodic orbits can simply be localized by 
extracting the positions of the fixed points from the TPM of higher iterates 
of the map. In our case of Fig.9 xpi and xpz are fixed points of / whereas 
{xf2, ^Fi} represents a period two orbit of /. 

The TPM of Fig.9 still contains much more information. Let us consider 
the vertical right and left edges of the individual objects which are part of 
the TPM. The positions of these edges are indicated in Fig.9 by yi, ?/2, y^. 
These quantities yi are the positions of the first or higher preimages of the 
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fixed point xp'?,- I/2, ?/3, 2/4 are first prcimagcs whereas 1/5,^/6 are second 
preimages of xpz- In general tliese vertical edges are given by the positions 
of the preimages of the fixed points with a positive derivative p'{xFi) > 
0. Approaching smaller and smaller values for x{n) we obtain increasingly 
higher preimages of the corresponding fixed points. After having discussed 
the positions of the individual objects in the TPM let us now turn to an 
investigation of the objects themselves. 

Apart from the objects of type (e) all constituents of the TPM in Fig. 9 
consist of a series of humps. The widths of the humps belonging to the objects 
of type (a) decrease strongly if the values of x{n) approach the corresponding 
values of the preimages iji of the fixed point xp-^. The objects (6) and (c) 
possess such an accumulation point of humps only on one side, namely at 
the position of the preimage ^2- In order to analyze this series of humps with 
decreasing widths we have illustrated in Figure 10 important parts of the 
objects labeled by (a), (6), (c), {d) in Fig.9. Figure 10 (a) shows the object 
(a) of Fig.9 in the neighbourhood of the preimage 1/3 of xp^^ i.e. 1) is 

illustrated as a function of — x{n) with a logarithmic abscissae. 

The eyecatching feature is the asymptotically (x(n) — > 1/3) exact scal- 
ing behaviour of the individual humps with a constant scaling factor. We 
therefore encounter not only the above-mentioned scaling behaviour of the 
series of objects of type (a), (e) but also a scaling behaviour of the interior of 
the objects themselves (apart from (e)). The corresponding interior scaling 
factor is in the case of Fig. 10 (a) P'{xp^), i.e. the derivative of the map 
P at the position of the fixed point {xpz) with a positive derivative. The 
lower turning points 5yyi for i = 1, 2, 3 indicated in Fig. 10 (a) represent a se- 
ries of increasingly higher preimages of the fixed point xp^ which possesses a 
negative derivative p' {xp^) < 0. This series converges in an asymptotically 
scaling way towards the preimage ys of 3:^3. The structure of the object (a) 
can be explained in even more detail. The series of humps can be obtained 
by intersecting increasingly higher iterates p"' of the map p with the line 
x{n -|- 1) = Xp4 and joining together the resulting humps. In this sense the 
object (a) describes the top x{n -|- 1) > Xi?4 of the humps of arbitrarily high 
iterates of the map p in the interval [y^, y^]. Most of the properties discussed 
for the example of the individual object labeled (a) can also be observed for 
the other constituents of the TPM in Fig.9. 

Figure 10 (b) shows the object labeled (b) in Fig.9 in the neighbourhood 
of the preimage 2/2, i-e. x{n -f 1) as a function of 1/2 — x{n) with a logarithmic 
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abscissae. The same asymptotic scaling behaviour now for x{n) y2 can be 
observed and the scahng factor is again determined by the stabihty of the 
fixed point {xpa) which possesses a positive derivative. The upper turning 
points dY2i for i = 1,2,3 indicated in Fig. 10 (b) represent a series of in- 
creasingly higher preimages of the fixed point which possesses a negative 
derivative P'{xp2) < 0. This series converges in an asymptotically scaling 
way towards the preimage y2 of 2:^3. The series of reversed humps can be 
obtained by intersecting increasingly higher iterates of the map p with the 
line x{n + 1) = xf2 and taking the lower part x{n + 1) < xf2- Figures 10 (c) 
and (d) illustrate the structures of the objects (c) and {d) in Fig. 9. Asymp- 
totic scaling properties are again observed and described by the stability of 
the fixed point with a positive derivative. For Fig. 10 (c) the lower turning 
points 5y2i for i — 1,2,3 which represent increasingly higher preimages of 
the fixed point Xp^ converge towards the preimage 1/2 of the fixed point Xps- 
The humps of the object (c) can be explained analogously to the humps in 
Fig. 10(a). Finally Fig. 10 (d) shows a scaling of reversed humps towards the 
preimage yi and another series of the preimages Syu of the fixed point xps 
which converges to yi. 

Let us briefly summarize the observations. Turning point maps belonging 
to higher iterates of the original map / contain a number of dynamically 
most relevant informations. It is easy to extract the position of not too long 
periodic orbits, their stability coefficients as well as the (higher) preimages 
of the fixed points of / by considering the TPM of the corresponding iterate 
The TPMs themselves as well as their disconnected constituents show 
an asymptotic scaling behaviour and exactly describe important parts of 
arbitrarily high iterates of the map /. 

To conclude our discussion of the TPMs of higher iterates /" we show 
in Fig. 11 the density Ptp2{x) of turning points for the doubly iterated map 
p. We observe step hke structures at the positions of the fixed points xpi 
which refiect the above-discussed properties of the corresponding turning 
point map. The singular behaviour of pTP2{x) at the boundaries of the unit 
interval is apart from a constant factor the same as for the density of the 
map itself. 
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5 Statistical properties of the turning point 
dynamics 

Having studied in detail the dynamical properties of the turning points in the 
previous chapters we turn now to an investigation of the statistical turning 
point properties. In contrast to the previous Sees, the results of the present 
Sec. hold only for doubly symmetric maps. A measure for the frequency of 
the turning points is given by the ratio of the number of turning points Ntp 
and the total number of points Ntoti i-e. Ptp = One half of that number 
(that is, only minima counted) has been called wave number [Ding, 1988] 
or over rotation number [Blokh and Misiurewicz] in the literature. Since 
we are dealing with a fully chaotic and ergodic system the quantity Ptp is 
solely determined by the necessary and sufficient condition x{n — l)>xi for 
a turning point x{n) and the frequency with which a certain region of the 
interval is visited. This frequency is however given by the invariant density 
Pf{x) of the map /. We therefore obtain the following expression for the 
portion of turning points of the chaotic trajectories 

Ptp= f\f{x)dx (22) 

J XI 

where xi = {1 — xp) is the inverse image of the fixed point xp > ^ of the 
map /. The question now arises how strong the quantity Ptp varies if we 
choose different maps /. Since all doubly symmetric maps are smoothly 
conjugate to each other we want to study the behaviour of the quantity Ptp 
under conjugation. Let the conjugation u transform the map g into /, i.e. 
/ = uogou~^. The invariant measure of g is then given by /Uc,(x) = iif{u{x)). 
It follows that the density transforms according to Pg{x) — pf{u{x))^ and 
we therefore obtain 

Ptv- l\, pMdy (23) 

The lower integration limit u~^{\ — xp) can be calculated in the following 
way. Since xp — fixp) — uo g o u~^{xf) we obtain u~^{xf) — g{u~^{xF)) 
which means that u~^{xf) is a fixed point of g. In addition the conjugation 
u possesses the property u{x) = 1 — ^(1 — x) which is due to the fact that 
we are dealing with symmetric maps. The image of the fixed point Xp of f 
is therefore also a fixed point of g, i.e. yp = u~^{xp) and we arrive at the 
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final equation 

Ptp= t Pg{y)dy (24) 

with yi — {1 — uf). This proves the invariance of the statistics of turning 
points under conjugation. All doubly symmetric maps possess therefore the 
same ratio of turning points and non turning points and Ptp is therefore a 
universal statistical quantity. Calculating this integral for the specific case 
of the logistic map yields Ptp = | which simply means that statistically two 
out of three points of the trajectories are turning points. 

6 Application of the turning point concept to 
the analysis of one dimensional time series 

The above-discussed characteristics of the turning point dynamics in the 
neighbourhood of the unstable fixed points of the chaotic system can in fact 
be used to analyze chaotic time series. The main goal of such an analysis 
is to find the location of the fixed points and higher periodic orbits of the 
underlying dynamical system. 

Since we are interested in fully chaotic and ergodic systems we encounter 
only hyperbolic fixed points (HFP). We have however to distinguish between 
HFP with and without reflection. For HFP with reflection subsequent it- 
erations of the map in the vicinity of the fixed point oscillate around the 
fixed point whereas for HFP without refiection a turning point in the neigh- 
bourhood of the fixed point is always followed by a stretching phase away 
from the fixed point. We will describe below an algorithm which allows us 
to determine the positions of the fixed points (periodic orbits) of a chaotic 
time series using these properties. The most important feature of our turning 
point method for the location of the fixed points is the fact that the above- 
described properties are not restricted to the linear neighbourhood of the fixed 
point. Therefore the position of the fixed point can be approximately deter- 
mined even in the case when the linear neighbourhood of the fixed point is 
not visited by the finite trajectory, i.e. time series, at all. 

As a first step of our algorithm we select those points of a given time 
series {x{i)\i — 1,...,N} which are turning points, i.e. which obey Eq.(l). 
The next step distinguishes between HFP with and without refiection. To 
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determine the HFPs with reflection we select those points {xt{i)} of the 
turning point trajectory (resulting from the first step) which have as a next 
iteration point in the original time series also a turning point {xt{i + 1)}. 
We group these points according to the values of their center of oscillations 
defined as Xcoii) = ^{xt{i) + Xt{i + 1)) . Thereby points with approximately 
the same vahics for Xcoii) belong to the same group. Using the amplitudes 
of oscillations d{i) = \xt{i + 1) — Xt{i)\ we look subsequently within each 
group for the turning point pair {xt{i),Xt{i + 1)) which corresponds to the 
minimal amplitude of oscillations. The position of the corresponding HFP is 
then given by the corresponding center of oscillation xp = ^{xt{i)+Xt{i + l)). 
By using the minimal amplitude we get rid of the big oscillations which occur 
due to the effects of the boundaries of the finite interval. 

We emphasize that the minimal amplitude of oscillations needs not to 
be smaller than the size of the linear neighbourhood of the fixed point in 
order to obtain a good estimation of the position of the fixed point. Our 
method possesses therefore an advantage compared to methods relying on the 
properties of the linear regime around the fixed point or recurrence methods. 
This is due to the fact that the turning around the fixed point is not a strongly 
localized property, and in particular not restricted to the linear regime. 

For the HFP without refiection we collect those points of the turning 
point trajectory which have as subsequent iteration points in the original 
trajectory non turning points, i.e. points belonging to a stretching phase of 
motion. If the stretching phase is increasing the value of the coordinate, the 
corresponding points belong to the group {xf{i)} and in the opposite case to 
{x^{i)}. As a next step we look for the extrema defined by Xp = 'min{xf (i)} 
and Xp = max{xt (i)}. The position of the fixed point is finally given by 

Xp — 2("^F Xp^. 

In order to demonstrate the efficiency of our algorithm we have applied 
it to a time series of the second iterate of the logistic map P{x) — f{f{x)) 
with f{x) = rx{l — x) for r = 3.9. Using a times series of only 200 points we 
are able to determine the positions of the fixed points of the system with a 
relative accuracy of 0.004. For the positions of the two HFP with refiection 
we get the values xpi — 0.36052 and xp2 — 0.89876 compared to the exact 
values 0.35897... and 0.89744..., respectively. For the HFP without reflexion 
we obtain xps = 0.74421 compared to the exact value 0.74359. An analysis of 
the same time series with, for example, recurrence methods yields a relative 
accuracy of 0.05 for the position of the fixed points and is therefore much less 
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accurate. Our algorithm for the turning point analysis of time series yields 
reliable results even if the number of points of the time series is comparatively 
small. 

The above algorithm can also be applied to the case of a scalar experi- 
mental ID time series, i.e. a chaotic signal, which originates from a higher 
dimensional dynamics. It is therefore by no means restricted to a ID time 
series arising from a ID system. The oscillating behaviour around the un- 
stable fixed points described above represent a characteristic feature of each 
component of a multidimensional signal. Indeed the chaotic dynamics of 
a multidimensional system in the neighbourhood of a fixed point can be 
characterized by the fact that the deflection or turning of the different tra- 
jectories approaching the fixed point varies significantly and therefore turn- 
ings can occur with respect to any direction of space. As examples for two- 
dimensional fully chaotic systems we mention the Ikcda [Ikcda,1979] or Henon 
map [Henon, 1976] which exhibits oscillatory behaviour with small amplitudes 
around the position of the fixed point for both components x as well as y. 

7 Summary and Outlook 

We have investigated the ergodic and fully developed chaotic behaviour of 
SSH maps in some detail. Our starting point was the observation that 
the chaotic trajectories suggest the existence of certain nontrivial structures 
around the unstable fixed point. A closer look at the trajectories showed us 
that these structures have their origin in the distribution of oscillations in 
the unit interval. A comparison of the densities for the center of oscillations 
as well as their amplitudes with the corresponding quantities for a weighted 
random map clearly revealed that the dynamics of the oscillations is highly 
organized and far from being random. Characteristics of the amplitude of 
oscillations as a function of the center of oscillations have been discussed. 

The central idea of our investigation is the introduction of a turning 
point analysis as a tool for the characterization of the dynamics of ergodic 
and fully developed chaotic systems. On our way of a quantification of this 
idea we have derived a turning point criterion which gives us a condition for 
the point x{n — 1) so that x{n) becomes a turning point. With the help of 
this criterion we were able to understand the dynamics of turning points. 
Roughly speaking the trajectories possess two phases of motion: a stretching 
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phase in the neighbourhood of the unstable fixed point which consists apart 
from its important first point only of non turning points and an oscillating 
stretching phase around the second unstable fixed point xp which consists 
exclusively of turning points. 

As an important tool we have subsequently introduced the turning point 
map (TPM) which ilhistrates the mapping of successive turning points. This 
TPM exhibits a number of dynamically interesting properties. It contains 
details of the first hump of iterates P for all i and maps in particular ar- 
bitrarily small iiiter\'als in the neighbourhood of the origin onto the finite 
interval [xi?, 1]. The TPM obeys an asymptotic scaling law and converges 
to a repeating functional form. A construction scheme for the humps in the 
asymptotic scaling regime has been derived. Another important property of 
the TPM is the fact that it allows us to locate the fixed point xp bs well as 
all its preimages. The turning point dynamics shows an infinite correlation 
length of the corresponding correlation function, i.e. the turning point anal- 
ysis extracts the correlated part from the originally completely uncorrelated 
trajectory. The origin of this strong correlation is the oscillating dynamics 
around the fixed point xp- 

We have shown that it is possible to define a so-called renormalized limit 
function which corresponds to the renormalized version of the infinitely iter- 
ated map. The renormalization is performed with respect to successive zeros. 
Successive humps are scaled onto the unit interval. We have illustrated the 
construction scheme of the limit function for the example of the logistic map. 

As a next step we have investigated the TPMs belonging to higher iterates 
of the original map /. The TPMs themselves as well as their disconnected 
constituents show an asymptotic scaling behaviour and exactly describe im- 
portant parts of arbitrarily high iterates of the map /. It is possible to 
extract the position of not too long periodic orbits of / by considering the 
TPM of the corresponding iterate A limited amount of information on 
the stability of the periodic orbits can be obtained from the scaling laws of 
the constituents of the TPMs. Informations on the position of high preimages 
of the fixed points can easily be extracted. 

We have performed a statistical analysis of the turning point dynamics 
for the special case of doubly symmetric maps. We could show that the ratio 
of turning points to the total number of points is an invariant for these maps 
and takes on the value |. This suggests that the ratio might be used for a 
characterization of nonsymmetric maps and divides the maps into different 
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equivalence classes. 

Our investigation and discussion of the dynamical and statistical be- 
haviour of the turning points, the TPMs as well as the correlation functions 
clearly show that a turning point analysis extracts valuable information from 
a given ergodic and fully chaotic system and is very well-suited for the anal- 
ysis of chaotic time scries. Using the main features of the dynamics of the 
turning points in the neighbourhood of the fixed point we developed an algo- 
rithm which allows us to locate the unstable fixed points and higher periodic 
orbits thereby distinguishing between hyperbolic fixed points with and with- 
out reflection. The most important property of our turning point method for 
the location of the fixed points is the fact that it is not restricted to the linear 
neighbourhood of the fixed point. Even in the case when the linear regime 
around the fixed point is not visited by the time series at all, it is possible to 
determine approximately the position of the fixed point. Our turning point 
analysis of the time series yields reliable results in particular if the number of 
points of the time series is comparatively small and represents therefore an 
attractive method for the investigation of dynamical systems of, for example, 
biological origin. 
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Figure 1: A chaotic and ergodic trajectory x{n) of the logistic map for 
1.5 X 10"^ iterations and the invariant density of the same trajectory for 
10^ iterations. 
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Figure 2: 50 iterations of a trajectory of the logistic map. The turning 
points are indicated with circles whereas the non-turning points are indicated 
with crosses. The stretching and oscillatory phases of motion are clearly 
visible. 
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Figure 3: The turning point density pTp of the logistic map. The step like 
structure is at the position of the fixed point xp. 
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Figure 4: (a) The density for the centers of the oscillations pLCo{x) of 
a typical chaotic trajectory of the logistic map. (b) The density for the 
amplitudes of the oscillations Plao{'^) of the same trajectory. 
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Figure 5: (a) The density for the centers of the oscillations pRco{x) of a 
trajectory of the weighted random map with Pr — Pl- (b) The density for 
the amplitudes of the oscillations Prao{,^^ of the same trajectory. 
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Figure 6: (a) The amplitudes of the oscihations as a function of the centers 
of the oscillations Xrao{^rco) foi' the random map. (b) The amplitudes of 
the oscillations as a function of the centers of the oscillations xlao{xlco) for 
the logistic map. 
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Figure 7: The turning point map of the logistic map for a logarithmic 
abscissae. Xj are successive preimages of the unstable fixed point a;^ = | on 
the left branch of /. 
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Figure 8: The absolute value |C(m)| of the correlation function C{m) — 
{{x — — {x)f))f for the logistic map. Circles indicate the corre- 

lation function of the turning points and crosses the correlation function of 
the original trajectory. 
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Figure 9: The turning point map of the doubly iterated logistic map for a 
logarithmic abscissae. The disconnected objects (a), (e) repeat themselves in 
an asymptotically scaling way if we approach the unstable fixed point 0. All 
objects show also an interior scaling behaviour (see Figure 9). yi correspond 
to the positions of the preimages of the fixed point with a positive derivative, 
i.e. preimages of xp^i = |. The positions of the fixed points xpi of p are 
indicated on the vertical axis. 
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Figure 10: Magnified parts of tlie objects (a), (b), (c), (d) of Figure 8 are 
shown on a logarithmic scale. The scaling behaviour of the interior of these 
objects becomes evident, (a) Scaling behaviour towards the preimage of 
xpshy a. series Sy^, dy32, (^2/33, •••of preimages of the fixed point xf4- (b) Scal- 
ing behaviour towards the preimage y2 of a;i?3 by a series 5l2i, ^^2, SY23, ... of 
preimages of the fixed point xf2- (c) Scaling behaviour towards the preimage 
y2 oixF3 by a series Sy2i, Sy22, Sy23, ... of preimages of the fixed point Xp4. (d) 
Scaling behaviour towards the preimage yi ofxps by a series Syn, dyu, Syi3, ... 
of preimages of the fixed point xf2- 
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Figure 11: The density ptp2 of the turning points of the doubly iterated 
logistic map p. Three step like structure are clearly visible. The positions 
of these steps are the fixed points xpi of p. 
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